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We present updated measurements of CP- violating asymmetries in the decays B° — > D* ± D^ 1 and 
B° -> D+D~ using (383±4) x W 6 BB pairs collected by the BABAR detector at the PEP-II B factory. 
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We determine the time-integrated CP asymmetry A d *± d =f = 0.12 ± 0.06 ± 0.02, and the time- 
dependent asymmetry parameters to be C D *+ D - = 0.18±0.15±0.04, S D ,+ D - — — 0.79±0.21±0.06, 
C D ,- D + = 0.23 ± 0.15 ± 0.04, S D ,- D+ = -0.44 ± 0.22 ± 0.06, C D+J3 - = 0.11 ± 0.22 ± 0.07, and 
S D + D - = — 0.54±0.34±0.06, where the first uncertainty is statistical and the second is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model (SM), CP violation arises from 
a complex phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark- mixing matrix, V Measurements of 
CP asymmetries in B° — > (cc)K^° decays Q by the 
BABAR 0] and Belle Q collaborations have firmly estab- 
lished this effect and precisely determined the parame- 
ter sin2,3, where (3 is arg[— KdV^/VtdV!^]. Another way 
to measure sin2/3 is to use decays whose amplitudes are 
dominated by a tree-level, color-allowed b — > ccd tran- 
sition, such as B° — > D ( -*^ ± D T . Within the frame- 
work of the SM, the time-dependent CP-asymmetries of 
B° — > D^^D^ are directly related to sin2/3 when cor- 
rections due to penguin diagram contributions are ne- 
glected. The penguin-induced corrections have been es- 
timated in models based on the factorization approxi- 
mation and heavy quark symmetry and are predicted to 
be a few percent [a, [||. However, contributions from 
non-SM processes may lead to a large shift 0. A signif- 
icant deviation in the sin2/3 measurement from that of 
the B° — > (cc)K^° decays would be evidence involving 
new physics beyond the SM. 

Studies of the CP violation in b — > ccd transitions have 
been carried out by both the BABAR and Belle collabo- 
rations. Most recently, the Belle collaboration reported 
evidence of large direct CP violation in B° — > D + D~ 
where C D+D - = -0.91 ± 0.23 ± 0.06 8], in contradic- 
tion to the SM expectation. However, such a large direct 
CP violation has not been observed in previous measure- 
ments with B° — » D^^D^T decays, involving the same 
quark- level weak decay 0,0, 11. 12|. 

In this Letter, we present an updated measurement of 
CP- violating asymmetries in the decays B° — > D* + D~, 
B° -> D*-D+ and B° D+D~. The datajised in this 
analysis comprise (383±4) x 10 6 T(45) -> BB decays col- 
lected by the BABAR detector at the PEP-II storage rings. 
The BABAR detector is described in detail elsewhere 1 1 31) . 
Monte Carlo (MC) simulation based on GEANT4 (lljis 
used to validate the analysis procedure and to study the 
relevant backgrounds. 

The decay rate /+(/_) for a neutral B meson decay 
to a common final state accompanied by a B a (B°) tag is 
given by 

f±(At) =e" |At|/r « /4r i3 o {(1 =p Aw) ± (1 - 2w)x 

[Ssm(Am d At) - C cos(Am d At)]} , (1) 

where At = t rcc — i ta g is the difference between the proper 
decay time of the reconstructed B meson (B rcc ) and that 
of the tagging B meson (B tag ), Tgo is the B° lifetime, and 
Arrid is the difference between the heavy and light mass 



eigenstates determined from the B°-B° oscillation fre- 
quency [IH . The average mistag probability w describes 
the effect of incorrect tags, and Aw is the difference be- 
tween the mistag probabilities for B° and B°. Since 
£)*+£)- and D*~D + are not CP-eigenstates, we can de- 
fine a time-integrated asymmetry Ad*±dt between the 
rate of B° -> D*+D~ and B° -> D*~ D+ , calculated as: 



A 



N 



D'±DT 



D* + D 



- — Nr>*- 



D*~D+ 



(2) 



N D , +D - +N D *- D + 

where N is the signal event yield. 

For B° — > D* ± D T , the general relations are Sd*±dt — 

— yjl — C 2 D , ±D ^ sin(2/3 c ff ± S), where 5 is the strong 

phase difference between B° — > D* + D~ and B° — > 
£)*—£)+ Under the assumption of negligible pen- 

guin contribution, f3 e g = (3, Ad*±e>^ = and C D ,+ D - = 
—Cjj*-u+. For B° — > D + D~ and in the case of negli- 
gible penguin contribution, C £>+£>- measures direct CP 
violation and is zero, while Sjj+jj- is — sin2/5. 

The selections of B° -> D*±DT and B° -> D+D~ can- 
didates are similar to those of our previous analysis [icj |. 
We reconstruct D* + in its decay to D°ir + . We recon- 
struct candidates for D° and D + mesons in the modes 
D° -> K-n+, K-n+ir , K-ir+ir+ir-, K°n + Tr~ 
D+ -> K-TT+TT+, K%+. We reconstruct B° 



and 
D+D- 

candidates only through the decay — » X =F 7r ± 7r ± . We 
require the reconstructed masses of the D° and D + candi- 
dates to be within 20 MeV/ c 2 of their respective nominal 
masses [TBI, except for the D a — > K~tt + it candidate, 
where we use a looser requirement of 40 MeV/c 2 . We 
apply a mass-constrained fit to the selected D° and D + 
candidates and combine D° candidates with a 7r + track, 
with momentum below 450 MeV/c in the T(45) frame, 
to form D* + candidates. 

We reconstruct the K° candidates from two oppositely 
charged tracks with an invariant mass within 20 MeV/c 2 
of the nominal K® mass [HI . The x 2 probability of the 
track vertex fit must be greater than 0.1 %. We require 
charged kaon candidates to be identified as such using a 
likelihood technique based on the Cherenkov angle mea- 
sured by the Cherenkov detector and the ionization en- 
ergy loss measured by the charged-particle tracking sys- 
tems [13}. We form neutral pion candidates from two 
photons detected in the electromagnetic calorimeter [l3| , 
each with energy above 30 MeV. The invariant mass of 
the pair must be within 30 MeV/c 2 of the nominal 7T° 



mass 
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and we require their summed energy to be 
greater than 200 MeV. In addition, we further apply a 
mass-constrained fit to the ir° candidates. 
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To suppress the e + e~ — > qq (q = u,d,s, and c) con- 
tinuum background, we exploit the contrast between the 
spherical shape of BB events and the more jet-like na- 
ture of continuum events. We require the ratio of the 
second to the zeroth order Fox- Wolfram moments [17J to 
be less than 0.6. We also use a Fisher discriminant, con- 
structed as an optimized linear combination of 11 event 
shape variables [181 ]: the momentum flow in nine con- 
centric cones around the thrust axis of the reconstructed 
B° candidate, the angle between that thrust axis and 
the beam axis, and the angle between the line-of-flight 
of the B° candidate and the beam axis. In addition, we 
employ a combined D flight-length significance variable, 
derived from the sum of flight lengths of the two D can- 
didates to reduce background. 

For each B° -» D^D* candidate, we construct a 
likelihood function £ maS s from the masses and mass un- 
certainties of the D and D* candidates [l9[ . The D mass 
resolution is modeled by a Gaussian whose variance is 
determined on a candidate-by-candidate basis from its 
mass uncertainty before the mass-constrained fit. The 
D*-D mass difference resolution is modeled by the sum 
of two Gaussian distributions whose parameters are de- 



ogous to that used in previous BABAR measurements de- 



termined from simulated events. The values of £„ 



and 



date energy Eg and the beam energy £Beam in the T(45) 
frame, are used to reduce the combinatoric background. 
From the simulated events, we optimize the maximum 
allowed values of — ln£ mass and |AJ5| for each individual 
final state to obtain the highest expected signal signifi- 
cance. 

We extract the signal yield from the events satisfying 
the selection criteria using the energy-substituted mass, 

toes = v^Beain — V*b i where P*b ^ s tne B° candidate mo- 
mentum in the T(45) frame. We select the B° candidates 
that have toes > 5-23 GeV/c 2 . On average, we have 1.5 
and 1.1 B° candidates per event for B° — > D* ± _D T and 
B° — > D + D~ respectively. If more than one candidate is 
reconstructed in an event, we select the candidate with 
the smallest value of — ln£ mass . Studies using MC sam- 
ples show that this procedure results in the selection of 
the correct B° candidate more than 95 % of the time. 

We perform an unbinned maximum likelihood fit to the 
Toes and At distributions to extract the CP asymmetries. 
We fit the events from B° -> D*+ D~ and B° -> D*~ D+ 
decays simultaneously. The probability density function 
(PDF) of the toes distribution consists of a Gaussian for 
the signal and a threshold function for the combina- 
torial background. We expect some background events 
to peak in the toes signal region due to cross feed from 
other decay modes. We estimate the fraction of events 
in the signal Gaussian due to this peaking background to 
be (8.8 ± 4.4) % for B° -> D^D* and (4.8 ± 7.4) % for 
B° — » D + D~ using detailed MC simulations of inclusive 
B decays. 

The technique used to fit the At distribution is anal- 



scribed in Ref. [2l|, |22j . We use information from the 
other B meson in the event to tag the flavor of the fully 
reconstructed B° -> D^D* candidate [Hj]. The sig- 
nal At PDF in Eq. [T] is convolved with an empirical At 
resolution function [21(. The At is calculated from the 
measured separation Az between the decay vertices of 
B 1CC and B tag along the collision (z) axis [2l|. The B tElg 
decay vertex is determined by fitting charged tracks not 
belonging to the B vec candidate to a common vertex, em- 
ploying constraints from the beam spot location and the 
-B rec momentum [2l[ . Only events with a At uncertainty 
less than 2.5 ps and a measured |Af| less than 20 ps are 
accepted for the fit to the At distribution. Both the sig- 
nal mistag probability and the At resolution function are 
determined from a large sample of neutral B decays to 
flavor eigenstates, B$ av . The combinatoric background 
At distributions are parameterized with an empirical de- 
scription that includes zero and non-zero lifetime compo- 
nents [2l|. The non-zero lifetime background is allowed 
to have effective CP asymmetries, and these float in the 
likelihood fit. By default, we assume that the peaking 
backgrounds have the same At PDF as the signal but 
zero CP asymmetries. 

The fits to the data yield 280 ± 19 signal events for 
B° -> D*+D-, 219 ±18 signal events for B° -> D*~ D+ , 
and 131 ± 14 signal events for B° — > D + D~, where the 
quoted uncertainties arc statistical only. In the region of 
™es > 5.27 GcV/c 2 , the signal purity is approximately 
41% for B° -> D*+D-, 34% for B° -> D*~ D+ , and 
46% for B° -> D+D~. The fitted CP violating parame- 
ters are 





0.12 ±0.06 ±0.02 


Cd'+d- = 


0.18 ±0.15 ±0.04 


Sd'+d- = 


-0.79 ±0.21 ±0.06 


Cd'-d+ — 


0.23 ±0.15 ±0.04 


Sd*-d+ = 


-0.44 ±0.22 ±0.06 


Cd+d- = 


0.11 ±0.22 ±0.07 


Sd+d- = 


-0.54 ±0.34 ±0.06, 



(3) 



where the first uncertainty is statistical and the second 
is systematic. 

Projections of the fits onto toes for the three different 
samples are shown in Figure [TJ Figure [2] shows the At 
distributions and asymmetries in yields between events 
with B° and B° tags, overlaid with the projection of the 
likelihood fit result. As a cross check, we repeat the fit by 
allowing the B° lifetime to float. The obtained lifetime 
is in good agreement with its world average [15 1. 

The systematic uncertainty of the time-integrated CP- 
asymmetry Ad*±dt is dominated by the potential differ- 
ences in the reconstruction efficiencies of the positively 
and negatively charged tracks (0.014). Other sources 
that contribute to the systematic error include the es- 
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FIG. 1: Measured distribution of m ES for (a) B° -> D*+ D~ , (b) B° -» D*~ D+ and (c) _B° -> candidates. The solid 

line is the projection of the fit result and the dotted line represents the background components. 



timate of the peaking background fraction (< 0.001), the 
uncertainty in the tties resolution for the B° — > D^D^ 1 
signal events (0.005), and a possible fit bias (0.004). 

The systematic uncertainties on C and S are evaluated 
separately for each of the decay modes. Their sources and 
estimates are summarized in Table Q] The systematic un- 
certainties arise from the amount of possible background 
that tends to peak under the signal and its CP asymme- 
try, the assumed parameterization of the At resolution 
function, the possible differences between the Bfl av and 
signal mistag fractions, the knowledge of the event-by- 
event beam-spot position, the uncertainties from the fi- 
nite MC sample used, the possible interference between 
the suppressed b — > ucd and the favored b — > cud ampli- 
tudes in some tag-side decays [23| , and the uncertainty in 
the Toes resolution for the signal events. All of the sys- 
tematic uncertainties are found to be much smaller than 
the statistical uncertainties. 

Since D* + D~ and D*~ D + are not CP-eigenstates, it 
is also illustrative to express the measured CP- violating 
parameters C and S in a slightly different parametriza- 
tion Q: C D > D = (C D , +D - +C D ,- D+ )/2, AC D * D = 
(Cd'+d- — Co*- d+) /2, Sd*d — {Sd'+d- + Sd*-d+)/2 
and ASd*d = {Sd*+d- ~ Sx>»- d+) /2. The quantities 
Cd* d and So* d parametrize flavor-dependent direct CP 
violation, and mixing-induced CP violation related to 
the angle (3, respectively. The parameters ACd'D and 
ASd* d are insensitive to CP violation. ACd*d describes 
the asymmetry between the rates T(B° — > D* + D~) + 
T(B a -> D*-D+) and T(B° -> D*~ D+) + T(B° 
D* + D~), while ASd*d is related to the strong phase 
difference, 6. We find 

C D , D = 0.21 ±0.11 ±0.03 
S D . D = -0.62 ±0.15 ±0.04 
ACd-d = -0.02 ± 0.11 ± 0.03 
AS D . D = -0.17 ±0.15 ±0.04, (4) 

where the first uncertainty is statistical and the second 
is systematic. 



In summary, this letter reports updated measurements 
of the CP violating asymmetries for the decays B° — > 
£)*±£)=f anc [ go £)+£)-. These measurements super- 
sede the previous BaBAR results [l(J, with a more than 
50 % reduction in the statistical uncertainties. The time- 
dependent asymmetries are consistent with the SM pre- 
dictions within their statistical uncertainties. We do not 
see evidence of large direct CP violation in the decay 
B° -> D+D- as reported by the Belle Collaboration @|. 
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0.026 


0.041 
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